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^ ; Abstract 

In nonrelativistic quantum mechanics and in relativistic quantum field theory, time t is a pa- 
rameter and thus the time-reversal operator T does not actually reverse the sign of t. However, 

O '■ 

l/^ . in relativistic quantum mechanics the time coordinate t and the space coordinates x are treated 



on an equal footing and all are operators. In this paper it is shown how to extend VT symme- 
try from nonrelativistic to relativistic quantum mechanics by implementing time reversal as an 
operation that changes the sign of the time coordinate operator t. Some illustrative relativistic 
►v>( I quantum-mechanical models are constructed whose associated Hamiltonians are non-Hermitian but 

C^ ' VT symmetric, and it is shown that for each such Hamiltonian the energy eigenvalues are all real. 



I. INTRODUCTION 

In nonrelativistic quantum mechanics the position x(t) is taken to be an operator while 
the time t is only a c-number parameter. To make quantum mechanics relativistic one must 
treat time and space equivalently. There are then two possibilities: one can either demote 
the spatial coordinates to parameters or promote the time coordinate to an operator. The 
former prescription is used in quantum field theory, where the field operators are treated as 
functions of the spacetime parameters x and t, but one can also construct sensible quantum- 
mechanical theories via the latter approach [l|, [j]- In such theories a new parameter is 
needed to parameterize evolution, and thus one introduces a fifth coordinate r that is an 
S0{3, 1) Lorentz scalar. In this five- dimensional formalism the space and time coordinates 
x^{t) become operator functions of r and one obtains an 5*0(3, l)-invariant relativistic 
first-quantized generalization of the nonrelativistic Heisenberg algebra [xj,^^] = iSj±: 

[x'^{t),p^{t)] = ^r/'^^ [x'^(r),p.(r)] = zS^^, (1) 

where t]^'^ is the 5*0(3, 1) Minkowski metric. 

The dynamics in this formalism is 50(3, 1) invariant in the four operators x^, but is 
nonrelativistic in the fifth coordinate r [because the dynamics is not 5*0(4, 1) or 5*0(3,2) 
invariant], and propagation is forward in r. However, just as the nonrelativistic quantum- 
mechanical operator x(t) can propagate forward and backward with respect to its time 
parameter t, in relativistic quantum mechanics all four components of x'^(t) can propagate 
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forward and backward in r [3]. The five-dimensional formalism of [1|, |2[] readily incorpo- 
rates forward and backward time propagation, so one can introduce antiparticles with first 
quantization alone without requiring the second-quantization techniques of quantum field 
theory. 

When the five- dimensional Hamiltonian operator H is Hermitian and its eigenfunctions 
have the separable form '?/'„(a;^,r) = (/)„(a;^)e~*'^"'^ and when its states obey the standard 
Dirac completeness relation 

J2\n){n\ = I, (2) 

the five-space forward propagator takes the form 

G5(x^,r;xr,0) = -t9ir){x'}\e-^''-\x':) = -z^(r) J] 0„(4)C(xr)e-^^-. (3) 



This propagator obeys a Schrodinger equation that is first order in r: 

[tdr + h) G^{x\ t- 0, 0) = 5{t)5\x^). (4) 

The four-dimensional propagators in the five- dimensional formalism are constructed by 
integrating out the fifth coordinate. Given ([3]), the associated four-dimensional propagator 
is then defined as 

/oo 
rfrG5(a;^,r;xr,0), (5) 

-CxD 

where TV is a normalization constant. Using the integral representation 



0{r) = -— du^, (6) 

one then obtains 

Ga{x'};x^) = N}^ _f I ,v • in 

Finally, since the wave functions are eigenfunctions of H, G^i^x'^] xf ) obeys 

-HG^{xf',0) = N6\x^'). (8) 

The primary objective in using the approach of l|, l2| is to choose a five- dimensional H 
so that Gi{x'^, 0) obeys a differential wave equation of the form 

D4G4{x^,0) = 6\x^), (9) 

where D4 is one of the familiar wave operators that appear in quantum field theory (such 
wave operators are typically higher than first derivative in time). [Normalizing (74 (x^, 0) ac- 
cording to ([9]) would fix the constant A^.] Thus, in the simple case where the five-dimensional 
Hamiltonian has the form H = p — {p^Y + ^1 © is the Fourier transform of the standard 
four- dimensional scalar field Feynman propagator D4 = d^d^. Because forward propagation 
in r gives the correct ie prescription for the causal Feynman contour in four dimensions, 
-D4 = dfj^d^ is the usual four- dimensional Klein-Gordon operator. Using the five- dimensional 
formalism, one can solve for a one-body quantum-mechanical Schrodinger-type propagator 
in five space, and from it one can construct a many-body quantum-field-theoretic propagator 
in four space. The five-space formalism also permits one to choose five-space Hamiltonians 
for which the resulting four-space propagator does not obey an equation of the form ([9]) with 
a familiar D4. In this paper we construct some simple models that lead to a propagator 
equation with a familiar 1)4 and some that have a more general structure. 



When Lorentz invariance was introduced in classical mechanics, it described the invari- 
ance properties of the line element ds^ = dt^ — dy?. In addition to invariance under the 
continuous orthochronous Lorentz transformations, the line element also possesses a set of 
discrete invariances, namely, space reflection P: x — ?► — x, t — ?► t, time reversal T: x — )■ x, 
t — 7- — t, and their product spacetime reflection VT: x — t- — x, t — ^ —t. However, when 
time reversal was introduced into quantum mechanics by Wigner, the time reflection of 
the id/dt operator was achieved not by replacing t by —t but rather by taking T to be 
an antiunitary operator that transforms i into ~i {T : i —^ —i); time t was treated as a 
c-number parameter that is not affected by T. In relativistic quantum field theory, time 
reversal is also not implemented by making the direct replacement t — )■ — t even though the 
line element ds'^ = dt^ — dy? possesses this time-reversal invariance. As noted above, in 
the five-dimensional relativistic quantum-mechanical approach used here, we treat time as 
an operator and thus we can implement a time-reversal operation that acts directly on the 
time. We can also implement VT transformations directly on the time operator. 

The ability to implement a VT transformation on the time operator is appealing because 
of the implications of VT invariance for Hamiltonians that are not Hermitian. In the last 
few years it has been recognized [^-13] that a quantum-mechanical Hamiltonian that is not 
Hermitian may still have an entirely real set of energy eigenvalues. In the cases that were 



explicitly considered in 



QJ? 



i-[3], the reality of the eigenvalues was traced to the existence of an 
underlying invariance of the Hamiltonian with respect to a combined VT reflection. Thus, 
while Dirac Hermiticity of the Hamiltonian is sufficient for reality of eigenvalues, it is not 
necessary. (Of course, Hamiltonians with entirely real eigenvalues can be both VT invariant 
and Dirac Hermitian.) However, recently it has been shown |8| that a Hamiltonian that 
is not VT invariant cannot have an entirely real set of energy eigenvalues. This means 
that VT invariance. in contrast to D.ao Herm,tic,ty. ,s necessary for the reality of energy 
eigenvalues [9| . (If one knows only that a Hamiltonian is not Dirac Hermitian, one can say 
nothing about the reality of the eigenvalues.) Thus, VT invariance of a Hamiltonian is a 
broader requirement than Dirac Hermiticity. 

In the non-Hermitian PT-invariant context we apply the five- dimensional formalism de- 
scribed above. To do this we recall [3] that when a Hamiltonian is VT invariant, its eigen- 
values are either real or they come in complex conjugate pairs. Consequently, both H and 
its Dirac-Hermitian conjugate W have the same eigenspectrum, and they are related by 



some similarity transform V 10 1 



VHV-^ = H\ (10) 

In this case, if \n) is a right eigenvector \R) of if, then {n\V rather then {n\ is a left eigen- 
vector (L| of H. Consequently, the energy-eigenstate-completeness relation ([2]) is replaced 
by 

Y.\R){L\ = T.W){n\V = I, (11) 

and ([3]) and ([7]) are replaced by 

G'5(x;;,r;xf,0) = -«^(r)(x^|e-*^"|xf) = -«^(r) ^(x^|n)e-'^""(n|V|xf ), (12) 

G'4(xy;xn = iV2^ F I • — • ^^^) 

— -C/„ -t- It 

The propagator (TT3|) is the relevant one in the VT case, and with its "PT-symmetric H, it 
also obeys ([8]). 

Invariance under VT reflection is a more physical requirement than Hermiticity because 
the proper orthochronous Lorentz group has a complex VT extension. Until now, this aspect 
of the Lorentz group has not been utilized because transformations that reverse the sign of 
the time have not been considered. In the present paper we consider such transformations 
and explicitly extend VT symmetry to the relativistic quantum-mechanical domain. In 
particular we study some simple non-Hermitian but PT-symmetric 5*0(3, l)-invariant model 
Hamiltonians using the five-dimensional formalism and for each Hamiltonian we show that 
all of the energy eigenvalues are real. 

II. A SIMPLE FIVE-DIMENSIONAL PT-SYMMETRIC HAMILTONIAN 

The generic five- dimensional action has the form / = /g^ (ir'L(r'), where r is the end 
point of integration. We begin with a simple example that illustrates the five-dimensional 
formalism. Specifically, we take a Lagrangian of the form 

2 
Tfl TTIUJ 

L = — x^i^ — {x^xf" - 2ia^x'' - a^a^) , (14) 

where fi = (0, 1, 2, 3), the dot denotes differentiation with respect to r, and a^ is a real, exter- 
nal, r-independent four- vector operator that commutes with x'^. As constructed, the action 



is a relativistic SO {3, 1) scalar function of the four x'^ coordinates, but it is nonrelativistic 
in the fifth coordinate r. We define a canonical momentum 

and then eliminate x'^ to obtain a canonical Hamiltonian 

H = p^x^ - L 

1 TflUJ 

= l^Pp^P^ + -Y' *^^^^^ ~ 2m^x'' - a^a^") . (16) 

The Hamiltonian flT6|) is not Dirac Hermitian because of the ia^x^ term. 

Next, we assign V and T quantum numbers to the x^ and p^ operators, just as we do 
with the nonrelativistic x and p = d'^/dt\ to wit, we take the three spatial components 
x^ to be V odd [Px'^(r)P~^ = —x'^{t)] and T even \Tx^{t)T~^ = a;'^(— r)], and take 
the three spatial components p'' = dx^/dr to be V odd [Vp^{t)V^^ = —P^{t)] and T 
odd [Tp^{t)T^^ = — P^(— t)]. Similarly, we take the time component x" to be V even 
[Vx^{t)V^^ = x^{t)] and T odd [T2;°(r)T^^ = — x°(— r)] and take the time component 
pO = rfxVc?^ to be V even [P/(r)p-i = /(r)] and T even [Tp°{t)T'^ = /(-r)]. With 
these assignments the four x^ are VT odd while the four p'^ are PT even. Because T also 
converts i to —i, these assignments are consistent with the commutation algebra in (Q. We 
summarize these assignments as follows: 



(17) 



Finally, we take the four-vector a^ to be VT even. For our purposes we will need a to 
be V even and thus T even, and a^ to be V odd and thus T odd. In the five-space the 
Hamiltonian f lT6|) is conjugate to r and not to x". Neither V nor T affect r because r is only 
a parameter, so with p^x^ = p^p^ /m being VT even, the Hamiltonian is VT symmetric. 

To determine the energy eigenvalues we take the spacetime metric to be diag(?7^;y) = 
(—1,1,1,1). Writing x^ = {t,x,y,z), we obtain a wave-mechanics representation of the 
algebra ([ID when p^ = —id/dx^; that is, 

■ d . d 

Po = -^^. V. = -^g^. (18) 
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p p*^ X x° 
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- + - + 


r 


- + + - 


VT 


+ + — 



Consequently, in five-space the Schrodinger equation takes the form 



rau 






The substitution y^ = x^ — ia^ brings (1191) to the form 



2m dxf^ dx'- 



iaf,){x 



^ - ian 



^{T,xn. 



(19) 



'' dr 



d' 



d^ 



+ 



niuo 



iy' - 1' 



Hr,yn, 



(20) 



2m \dt'^ 9y2 

and reduces the Schrodinger equation to a four-dimensional harmonic oscillator with 
Minkowski signature. Noting that 



d 



d 



rau 



T^^-^ +^(y -^^) 



-mu>{y^-t^)/2 ^ 2cue"™'^(y'~*'^/^ 



(21) 



2m yOf^ dy^J ' 2 

we see that the t-dependent sector contributes a positive zero-point energy equal to uj/2 just 
as the y-dependent sector does. Because all the eigenvalues of a harmonic oscillator are real, 
the five-space energy eigenvalues of (HM are given by 



E. 



[n. 



c + Uy + Hz + fit + 2)UJ, 



(22) 



where each of n^, Uy, n^ and n^ ranges over the positive integers. Thus, while the Hamiltonian 
( IT6l) is not Hermitian, all of its energy eigenvalues are real. 
For this model the five-space propagator obeys 



' d I d d 



rau 



dr 2m dx^ dx^ 



{Xf, -ia^){x^ -ia^) 



G'5(x^r;0,0) = 5(r)5^(x'^) (23) 



and we show in Appendix A that 



G5(x^r;0,0) = ^(r) 



sm ujT 



exp 



imuj cos ujt{x^ — ia^) {x^ — ia^) 



2sinajr 



(24) 



The propagator of the associated four-dimensional theory is then obtained via ([5]), and it 
obeys (|8]) with the "PT-symmetric H = —d^d^/2m + muj'^{x^ — ia^){x^ — ia^^)/2. 

Using the "PT-theory techniques described in |6| , one can demonstrate the reality of the 
eigenvalues algebraically without actually solving the Schrodinger equation. To do so, one 
must construct an operator e^ that possesses four key properties: (i) a similarity transfor- 
mation using e® preserves the commutation relations; (ii) Q is a Hermitian operator (so 
that e^ is not unitary); (iii) like V in flTO|) . e^ effects the transformation 



e-^He^ = H^; 
7 



(25) 



(iv) the operator 

H = e-2/2i7eS/2 (26) 

obeys W = H. The existence of such a Q operator imphes that the energy eigenvalues of 
H are all real. 

We now construct the Q operator for our simple five-dimensional model. Note that the 
momentum operator will effect the transformation 

e-^'^P^x^e^'P" = x^' + ib^", (27) 

and leave the commutation relations ([1]) untouched for any four-vector b'^ that commutes 
with both x^ and p^. Given fl27|) . we identify Q as the Hermitian operator 2a'^p^ because 



Similarly, the transformation 

1 7T)(jJ ~ 

e-^^P^He^'P" = p^pf" + x.x" = H (29) 

2m 2 

generates an equivalent Hamiltonian H that is manifestly Hermitian. 

In VT quantum mechanics one introduces an operator C that is required to obey 

[C, H] = 0, C^ = I. (30) 

One constructs this operator by making the ansatz C = e^V, where the operator V obeys 
V^ = I. In this form, the operator C fulfills the condition C^ = I provided that Q satisfies 
VQV = -Q. With e-2 generating e'^He^ = W, the operator C obeys C~^HC = H ii V 
generates VHV = H"^ . For the Q and H of interest here, both VQV = —Q and VHV = H^ 
hold provided that aP is V odd and a is P even. With this choice for the parity of a'^, we 
then identify C = e^V. (Previously, we had required that a^ be Vl~ even.) Then, if both a'^ 
aud p- a..e VTeven, the operator Q is VT everr. As eonstrueted. C thus obeys |C. PT] = 0, 
as expected [8|, |9[ when all energy eigenvalues are real |11| . 

III. FIVE-DIMENSIONAL PAIS-UHLENBECK OSCILLATOR 



In 1950 Pais and Uhlenbeck 12| explored the question of whether the Pauli-Villars reg- 



ulator associated with the fourth-order equation of motion 

(^2 _ V2 ^ ^2) (^2 _ ^2 ^ M|)0(x,t) = (31) 



and propagator 

^^^'^ " (A;2 + M2)(A;2 + M|) " M[^^ (p + Mf ~ PTm|J ' '^^^'' 

where /c^ = — (A;°)^ + k , could be physically viable, or whether it was merely a mathematical 
technique to regulate Feynman integrals. To this end they replaced the scalar field 0(x, t) 
by a single coordinate z{t) and examined single momentum modes w^ = k^ + M^ and 
u;| = k^ + M|. The resulting nonrelativistic quantum-mechanical limit of the equation of 
motion (13T1) and the propagator ( l32l) . 

(9?+.?)(a? + .S.W = 0, G(E) = -^(^-^), (33) 

is known as the PU oscillator. 

Pais and Uhlenbeck found that if the theory were quantized with a standard positive- 
metric Hilbert space, the energy spectrum would not be bounded below. One can evade this 
negative-energy problem by quantizing the theory in a negative-metric Hilbert space, but 
as the relative minus sign in f l33p indicates, the disadvantage of doing so is that one obtains 
states of negative Dirac norm and evidently loses unitarity. 

The PU oscillator was revisited in 2008 [l3|, [l^ and a new realization of the theory was 
found in which the Hilbert space has neither negative-energy nor negative-norm states. In 
this realization the Hamiltonian is not Dirac- Hermitian but is instead VT invariant. The 
norm is given by {L\R) = (?2|V^|n), rather than by the Dirac norm {n\n), and the complete- 
ness relation is given by ( ITTl) rather than by ([2]). In analogy with ( TT3l) . the relative minus 
signs in ( l32l) and ( l33l) are generated by the presence of the V operator in the propagator 
and not by quantizing with an indefinite metric. This realization took a long time (more 
than half a century) to discover because the Hamiltonian of the theory appeared to be Dirac 
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Hermitian even though it is not. (In Refs. [13|, |1J] the nonrelativistic Vl~ realization of the 



PU oscillator is studied, and in Ref. [IJ] the relativistic scalar field theory is examined.) 

For the case of the nonrelativistic PU oscillator, the equation of motion ( l33l) for the 
coordinate z{t) can be derived by a stationary variation of the PU oscillator action 



hv = ^ I dt 



•■2 12: 2\ -2 I 222 
Z —{UJi+Ul2]Z + UJ1UJ2Z 



(34) 



where 7, Ui and U2 are positive constants. Since z serves as the conjugate of both z and z, 
the action is constrained. One thus replaces i by a new variable x, and using the method 



(35) 



of Dirac constraints, one obtains 15|, ll6| the Hamiltonian 

2 

Hpu = ^ + p^x + - {ul + ujfj x^ - -^1^2^^ 

with two canonical pairs that obey [x,px] = i and [z,Pz] = i- 

In the reahzation of the theory for which the energy eigenvalues are bounded below, -ffpu 
appears to be Hermitian but it is not. Specifically, one solves the Schrodinger equation for 
the ground state of the system with energy Eq = {ui + W2)/2. The eigenfunction is 



ipo{z,x) = exp 



-{Ui + U2)uJlUJ2z'^ + i'yU)lUJ2ZX - -{ui + UJ2)x^ 



(36) 



This eigenfunction diverges exponentially for large z, so integration by parts generates sur- 
face terms that cannot be discarded. Thus, one cannot represent the operator p^ by —idz. 
However, one can replace z by iz (this is equivalent to working in a Stokes wedge in the 
complex- 2; plane that includes the imaginary z axis but not the real one ll3|), and represent 
Pz by —idiz = —dz- The eigenfunction then vanishes exponentially as z becomes large. The 
highly unusual implication of the structure of f p6|) (and the reason it took so long to find) is 
that while both conjugate pairs of coordinates are obtained from the same Lagrangian, the 
commutator [x,^^.] = i is realized by Hermitian operators, while the commutator [z,pz] = i 
is realized by anti- Hermitian operators. As a result, the pzX cross-term in (l35l) is not Her- 
mitian, and the Hamiltonian Hpjj is also not Hermitian. 

Rather than using non-Hermitian operators, we make the similarity transformation 



y 



^npzz/2^^-np2z/2 



-IZ, 



Q 



^■Kp^z/2 ^-■Kp^z/2 



iPz 



(37) 



to construct Hermitian operators y and q that obey [y,q] = i. In terms of y and q the 
Hamiltonian now takes the form 

Hpu 



P ■ , 7 / 2 , 
--^qx + -[u,+ 



2\ 2,7 222 

^2 j ^ + 2^1^22/ > 



(38) 



where for notational simplicity we have replaced p^ by p. The Hamiltonian H^u is now 
manifestly non-Hermitian. 

While H-pu is not Hermitian, the V and T quantum-number assignments 



(39) 





p X q y 


V 


+ + 


r 


- + + - 


vr 


+ - + - 
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make Hpu symmetric under VT reflection. Introducing the operator 

Q = apq + Pxy, a = log ( ) , /3 = a^y'^ufu'^, (40) 

'JU1U2 \UJi— UJ2/ 



m 



we then find that |13l . Il4| the similarity-transformed PU Hamiltonian 

represents two uncoupled harmonic oscillators. The transformed Hamiltonian -ffpu in ( HTl) 
is both Hermitian and manifestly positive definite. This realization of the quantum theory, 
which is associated with the non-Hermitian Hpu, has no negative- norm or negative-energy 
eigenstates [l7|. 

Because the transformation with e^/^ is not unitary, the propagator 

D{Hpy) = (x', y'le-'^'^^'lx, y) = (x', y'|e2/2g-^i/putg-Q/2|^^ y^ (42) 

associated with ifpu does not transform into the propagator 

D(^Pu) = (x',y'|e-^^™*|a;,2/) (43) 

that one would ordinarily associate with a two-uncoupled-oscillator system. The state 
{x,y\e^/'^ is not the conjugate of e~®/^|a;,y), and the propagators in f H2]) and fH51) are 
not equivalent; for this realization of the PU Hamiltonian we must use f l42p and not (H5|l . 
The dependence on the operator V^ = e~^ is crucial because it generates the relative minus 
sign in fl33|) . 

We now illustrate VT invariance in relativistic quantum mechanics by applying the five- 
dimensional formalism to the PU oscillator. We will see that a straightforward covariant 
generalization of the PU oscillator does not lead back to f l32p . Consequently, in the next 
section we provide an alternate five-dimensional formalism that does. 

To generalize the PU oscillator to relativistic quantum mechanics we replace (IM|) by 

1=1 fdr [z^z'' - (m2 + Mi) z^z^ + MlMlz^z^'] , (44) 

where the dot denotes differentiation with respect to r. Because of constraints associated 
with this action, the Hamiltonian has the form 

H = i^^)|M: + (p^)^^. + 1 (m,2 + Mi) x,x^ - lMlMlz,z>^. (45) 

11 



Recalling the transformation in ( 137|) . we let y^ = —iz^ and q^ = i{pz)^- On setting 
{Px)^ = P^ we obtain two canonical pairs of operators that obey 



[x'^(r),p^(r)] = ^r/'^^ [q>^ (r) , y^ (r)] = ii^^\ 



and a Hamiltonian of the form 



H 



27 



iq^x^ + I [Ml + M^) x^x^ + ^-MlMly.y^^. 



(46) 



(47) 



The assignments 





p p° X a;° q g° y y° 


V 


- + - + + - + - 


T 


- + + - + + 


VT 


+ + + + 



(48) 



in which x" changes sign under T, then establish that -ffpu is VT symmetric. 
Next, we introduce the operator 

1 



Q = apf.q'' + idXf.yf' 



a 



lo iMi±^ 



7M1M2 VM1-M2 



and find that 



H = e-Q/'He^/' = ^ 



Qf.q'' 



I Ti/r'i 



(3 = aYMfM^ 



' 7i/r2 7i/r2„ 



Mfx^x'' + -^MfM^y^y^^. 



(49) 



(50) 



27 ' 27Mf ' 2'"'^"'^'" ' 2' 

Thus, the energy eigenvalues of the PT-symmetric Hamiltonian H are all real. 

We show in Appendix A that if we set Mi = M and M2 = 0, the five- space propagator is 

gJB/A 



G5(x^l/^r;O,O,O) = 0(r) 



A2 



(51) 



where 



25/7 = Mx^x'' (sin Mr - Mr cos Mr) -M^y^y^ sin Mr + 2zM2x^y^(l- cos Mr), 
A = 2 - 2 cos Af r - Mr sin Mr. (52) 

The propagator of the associated four- dimensional theory may now be obtained by 
performing the integral in ([5]), and the resulting propagator will obey ([8]) with H = 
— (l/2'y)d/dx^d/dx^ — x^d/dy^ + 'jM'^x^x^/2. While of interest in itself, this propaga- 
tor is not of the generic Pauli-Villars form given in (132|) . Thus, in Sec. IIVI we provide an 
alternate choice for the five- dimensional Hamiltonian that will lead to ( 132|) . 



12 



IV. ALTERNATE FORMULATION OF THE FIVE-SPACE PU OSCILLATOR 

Given the structure of (!3T|) we take the five-space H to have the operator form 

H = -[-(p^+f + Mni-iP'f+f + M',]. (53) 

For this Hamiltonian the five-dimensional energies are given by 

E, = -[-{p'Y+p' + Ml][-{py+p' + Mi], (54) 

where the momenta in fl54p are the eigenvalues of the operators in (IS5]) . Inserting these 
energies into ([7]), we obtain the Pauli-Villars propagator in ( 15^ . with i^ being satisfied. 

Equation ( 15^ leads directly to ([22]), but its use here is nonstandard because it does 
not have a simple Lagrangian counterpart. In the previous examples and in the deriva- 
tion of the Klein-Gordon propagator, one can start with a five-dimensional action (of the 
form Jq drx^xf^ for the specific Klein-Gordon case) and by a canonical procedure derive 
a Hamiltonian from it. The Lagrangians in these examples are quadratic functions of 
the coordinates, so the procedure is straightforward and yields Hamiltonians that are also 
quadratic. However, the Hamiltonian f lS^ is not quadratic; it is quartic because the wave 
operator in fl3T]) is a fourth-order derivative operator [18]. Since the Lagrangian is given by 
L{x^) = p^x^ — H{p^p^) and since p^ = dL/dx^, one can in principle construct L{x^) if 
one knows H{jp^p^). Doing so for ( 153|) is difficult, so we start directly with H{p^p^). Once 
we have H{p^p^), we can then use the representation in ([71) without needing to know the 
structure of the Lagrangian. 

We can recover the four-dimensional Pauli-Villars propagator, but at first it appears that 
the Hamiltonian in (153|) is Hermitian. Moreover, in the second-order Klein-Gordon case 
with H = —{p^y + p^ + Mf and real E^ the Hamiltonian is Hermitian. However, in the 
fourth-order case, we note that (p°)^ is given as 

ipr = \ (El + El ± m - Elf - AE,YI') , (55) 

where Ef = p^ + Mf. Thus, now there can be real values of E^ for which (p°)^ is complex 
and for which the operator (p°)^, and thus H, is not Hermitian. (Note that with E^ being 
real, the Hamiltonian must be VT invariant.) 

In addition, we note that for general Mi and M2, if we take E^ to be zero, the eigen- 
functions associated with the operator H in ([53]) will have the form ipi = q-^^i^+w-^ and 
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'?/'2 = e~*'^^*^*^'^. However, if we then set Mi = and M2 = 0, there will be eigenfunc- 
tions of the form ipa = e"*^*"*"*^'^' and tpb = n^x^e"*^*"''*^'^, where n^ is the unit timelike vector 
n^ = (1, 0, 0, 0). Of the two ipa and ipb eigenf unctions, only ipa is stationary; ipi, grows linearly 
in the time coordinate, which indicates that the Hamiltonian has Jordan-block form and that 
it has an incomplete set of eigenvectors. Consequently, the Hamiltonian H in (1531) cannot 
be diagonalized and is not Hermitian. Since H is not Hermitian when E^ = Mi = M2 = 0, 
it must also not be Hermitian for a range of values of these parameters. In Ref. [ij] it was 
found that in the equal-frequency limit ui = 002 of the PU oscillator, the Hamiltonian in 
fl55]) is also nondiagonalizable and non-Hermitian. 

The solutions to f l55p thus break up into two sectors. In one sector the Hamiltonian 
is Hermitian and the energy eigenvalues are unbounded below {AE^ < (M^ — M'^Y). In 
the other sector the Hamiltonian is not Hermitian and the energies are bounded below 
(4£'5 > {Ml — M|)^), just as in the case of the nonrelativistic PU oscillator. If E^ is real, 
the Hamiltonian is VT invariant in both cases. In the sector where H is Hermitian the four- 
space propagator is given by (JTj). In the non-Hermitian sector the four-space propagator 
is given by ( TT3l) and as before, the V operator then generates the relative minus sign in 



the Pauli-Villars propagator 19| . Our five-space treatment of the Pauli- Villars propagator 
based on (l53l) recovers the key features of the analyses of Refs. [l3|, llJl- We see that one 
can extend VT symmetry to the five- dimensional formalism, and while we have not directly 
studied the time-reversal and VT properties of the time operator in the Pauli-Villars case, 
those properties follow directly from the commutation relations ([I]) depending on how they 
are explicitly specified for (p)°. 

V. SUMMARY 

Using a number of elementary models, we have shown in this paper that the standard 
techniques of VT quantum mechanics extend and apply to relativistic quantum mechanics, 
where the time- reversal operator T reverses the sign of the time operator x°. We conclude 
that relativistic PT-symmetric quantum mechanics is physically viable. 

The work of CMB is supported by a grant from the U.S. Department of Energy. 
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Appendix A: Construction of Five-space Propagators 

To construct propagators that obey the five-dimensional equation 
(idr + H) G5{x^,t;0,0) = 6{t)5^{x'^), we first recall how a propagator is constructed 
when the eigenmodes of H are plane waves. For the nonrelativistic quantum-mechanical 
free particle in one space dimension there is a plane wave basis and the propagator is given 
by 

Gi{x,t;0,0) = -^ /t/pe-'P"-'P'*/™. (Al) 

ZTT J 

When idt acts on —i9{t), we generate the 6{t)6{x) term, while if we omit the 9{t) function, 
the rest of the propagator obeys 



' d 1 &■ 



2 



dt 2m dx 



Ri{x,t;0,0)=0, (A2) 



where G'i(x, t;0,0) = 9(t)Ri{x,t;0,0). The Fourier transform in (1A1|) can be performed 
analytically and yields 

Gi(x,t;0,0) =0(t) (^£^y^%— V2^ (A3) 

The term /stat = rnx'^/2t in the exponent is the value of the classical action / = 
(m/2) JqcHx'^ for the stationary path x = between the end points (x = 0,t = 0) and 
{x,t). 

If one were to calculate this propagator as a path integral J[dx]e^^ over a complete basis 
of paths between the end points, one would obtain the same e*^^'^^'^ phase, but one would 
not know the multiplicative pre-factor. This pre-factor is determined by requiring that 
the propagator obey ( 1A2I) . (If one does not have a plane- wave basis, one can evaluate the 
propagator via a path integral and then use the Schrodinger equation to determine the 
pre-factor.) 

For the one-dimensional harmonic oscillator (where the basis is not plane waves), the path 
integral again has the form c^^^tat^ where /stat is the value of / = (m/2) /q (it[x^ — w^x^] as 
evaluated in the stationary path x + u^x = between the end points {x = 0,t = 0) and {x = 
Xf,t = T). Noting that x^ — u'^x'^ = d{xx)/dt — xx — u'^x'^, we obtain /stat = mxfXf/2. The 
solution to the equation of motion is x{t) = Xfsinut/ sincoT, x(t) = uxj cos out/ sinuT, so 
we obtain /stat = ma;x^cosc<;T/2sina;T. With this form for /stat, the pre-factor evaluates 
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to (sin uT) ^/^ and the propagator is 



^ , r^ X ^/^N / 1 \^^^ fimcox'^ COS ujT\ ,. ,, 

G^{x,T;0,0)=e{T)(——) exp . A4 

\smuTJ \ 2smuT I 



The propagator (HM is the shifted covariant generahzation of this result. 

The propagator associated with the PU oscillator action given in ( 1M|) has already been 
reported in the literature |20|, and because the action is quadratic, the J d[z] path integral 
between end points with fixed z and z has the form exp(iJsTAT) with the appropriate Jstat- 
Here, we present a simplified version of the propagator in which we set ui = u, U2 = 0. In 
this case the classical action reduces to 

Ip^ = ^Jdt{^-u'z'), (A5) 

and the stationary classical equation of motion is given by 

d^{z + coh) = 0. (A6) 

Noting that 

dt (zz - zdlz - uj'^zz) =z^ - u^z^ - zdf {z + uj'^z), ( A7) 

on evaluating Jstat between z = 0, z = at t = 0, and z{T), z(T) at t = T, we obtain 

/sTAT = (7/2) {ziTYziT) - z{T)d^z{T) - u'z{T)z{T)) . (A8) 

Hence, introducing 

coaA(T) = z(T){ijjT — sin uT) — (jjziT) {1 ~ cos uT), 
(3A{T) = z{T){l - cos uT) - z{T)uj sin uT, 
A{T) = 2-2cosLoT-ujTsinujT, (A9) 

we find that the solution to flA6l) that satisfies the boundary conditions takes the form 

z(t) = —a{l — cosijjt) — {f3/ijj)sinut + (3t, 

z{t) = —au sin out — (3 cos ut + f3, 

z(t) = —au'^ cos ut + (300 sin ut, 

dfz{t) = au^ sin ut + (3u^ cos ut. (AlO) 
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In this solution /stat obeys 

--^STAT = u)z'^{T) {sin uT — uT cosuT) — 2u''z{T)z{T){l — cosuT) + uj'^z^{T) sinwT. 



7 

(All) 

Finally, we verify that this function is a solution to the Schrodinger equation associated with 

( l35l) and identify the pre-factor as A~^^'^{T). The propagator is thus A~^^'^{T)e^^^'^^'^ . Its 

covariant generalization, obtained by using ( 1371) . is given in ( 15T1) . 



[1] R. P. Feynman, Phys. Rev. 80, 440 (1950). 
[2] Y. Nambu, Prog. Theor. Phys. 5, 82 (1950). 
[3] For a more recent discussion of propagation and path integration in the five-dimensional 

formalism see e.g. P. D. Mannheim, Phys. Lett. B 137, 385 (1984); Phys. Rev. D 32, 898 

(1985); and Phys. Lett. B 166, 191 (1986). 
[4] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 5243 (1998); C. M. Bender, S. Boettcher, 

and P. N. Meisinger, J. Math. Phys. 40, 2201 (1999). 
[5] C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. Lett. 89, 270401 (2002); 

Am. J. Phys. 71, 1095 (2003); Phys. Rev. Lett. 93, 251601 (2004). 
[6] C. M. Bender, Contemp. Phys. 46, 277 (2005) and Rep. Prog. Phys. 70, 947 (2007). 
[7] P. Dorey C. Dunning and R. Tateo, J. Phys. A 34, L391 (2001) and 34, 5679 (2001); 



Czech. J. Phys. 54, 35 (2004); |arXiv:hep-th/0201108t J. Phys. A: Math. Gen. 40, R205 
(2007). 

[8] C. M. Bender and P. D. Mannheim, Phys. Lett. A 374, 1616 (2010). 

[9] Technically, what was shown in 8|] is that the secular equation that determines energy eigenval- 
ues is real if and only if the Hamiltonian is VT symmetric. Since the solutions to an equation 
that is not real cannot all be real, VT invariance of a Hamiltonian is thus a necessary condition 
for the reality of eigenvalues. The solutions to a real equation can either be real or appear in 
complex conjugate pairs. Thus, VT invariance alone does not guarantee reality. However, as 
g], whether one has the real case or the complex-conjugate pair case depends 



noted in Ref. 



6|] commute with 



on whether or not the available C operators of VT theory described in Ref. 

VT- Commutation of both the Hamiltonian and all available C operators with VT is then 

both necessary and sufficient for the reality of energy eigenvalues. 

17 



[10] In P. D. Mannheim, PT symmetry as a necessary and sufficient condition for unitary time 

evolution, larXiv:0912.2635K ^l [hep-th] and A. Mostafazadeh, Int. J. Geom. Meth. Mod. Phys. 

7, 1191 (2010) the converse was shown: namely if there exists a V that effects VHV^^ = H^ , 

then H must be VT symmetric. 
[11] While we do not discuss it here, we note that analogous to the Hamiltonian given in (J16p . 

comparable results can be obtained for the Hamiltonian H = {l/2m){p^ — iCfj){p^ — ic^) + 

muP'x^x^^ /2, where c'^ is a real "PT-odd four-vector. 
[12] A. Pais and G. E. Uhlenbeck, Phys. Rev. 79, 145 (1950). 
[13] C. M. Bender and P. D. Mannheim, Phys. Rev. Lett. 100, 110402 (2008); J. Phys. A 41, 

304018 (2008). 
[14] C. M. Bender and P. D. Mannheim, Phys. Rev. D 78, 025022 (2008). 
[15] P. D. Mannheim and A. Davidson, |arXiv:hep-th/0001lT5| 
[16] P. D. Mannheim and A. Davidson, Phys. Rev. A 71, 042110 (2005). 
[17] Transforming -ffpu in ([55]) with Q = iapxPz — if3xz to H = p'^/2^ + 7a;fx^/2 — p^/27a;f — 

'yuJiUJ2z'^ /2, we see that if pz and z were Hermitian, the Hamiltonian would be unbounded 

below. 
[18] One can either represent the fourth-order derivative theory as a quadratic theory based on two 

sets of conjugate position and momentum pairs or as a quartic theory based on one conjugate 

pair alone. Both representations yield the same energy eigenvalues. 
[19] A relative minus sign can appear in (j32p in the Hermitian Hamiltonian case where there are 

no states with negative norm because as noted in \l4\ . in this case the contour for the p^ 

integration is not the standard Feynman one. Specifically, in the generic form 1/{A + ica) — 

1/{B + ieb) = {B — A)/(AB + ie), where e = Bea + Aef, the sense of the contours depends on 

the values of A and B. 
[20] P. D. Mannheim, Found. Phys. 37, 532 (2007). 



18 



